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ABSTRACT 

A linear relation between the hydrogen column density (Nh) and optical extinction 
(Ay) in the Galaxy has long been observed. A number of studies found differing results 
in the slope of this relation. Here, we utilize the data on 22 supernova remnants that 
have been observed with the latest generation X-ray observatories and for which optical 
extinction and/or reddening measurements have been performed and find NH(cm^^) = 
(2.21 ± 0.09) X 10^^ Ay (mag). We compare our result with the previous studies and 
assess any systematic uncertainties that may affect these results. 
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1 INTRODUCTION 

Photoelectric absorption by interstellar matter causes rapid 
attenuation of the observed soft X-ray spectra in all Galactic 
sources. Measuring the amount of the X-ray extinction yields 
information on the total column density of the interstellar 
matter between the observer and the source, because any 
element that is not fully ionized contributes to the extinc- 
tion of X-rays. Although the amount of the X-ray extinction 
is generally expressed in terms of the equivalent hydrogen 
column density (Nh), it is generally caused, especially above 
~0.25 keV, by the most abundant heavier elements like O, 
Ne, Fe, Mg, and Si. 

Optical extinction, on the other hand, is caused by 
■ grains composed of these same heavier elements. Therefore, 
the relation between optical and X-ray extinction depends 
on the relative abundances of metals and the temperature 
of the intervening medium, which primarily determine the 
relative fraction of metals in grains versus in the gas state. 
For sources, however, that lie at large distances from the 
observer, photons traverse many different regions of the in- 
terstellar medium (ISM), and, thus, sample a wide variety 
of physical conditions. Under the assumption that different 
Unes of sight sample the same distribution of physical condi- 
tions in the ISM, one expects a relation between the optical 
extinction and the hydrogen column density towards distant 
sources. 

The relation between optical extinction and hydrogen 
column density has been observationally studied for decades 
using various te c hniqu es (see, e.g., Reina fc TarenghH 
1 19731 : iGorensteinI Il975l : iPredehl fc Schm itt 199^). Us- 
ing measurements for two X-ray binaries together with 
the two ext ended sources GCX a nd Cas A (SNR 
G111.7-02.11. lReina fc Tarenghil (j 19731 ) derived the first re- 
lation NH(cm"^) = 1.85x 10^^ Av(mag). [Gorensteinl (119751 ) 



found a similar relation shortly afterwards using indepen- 
dent optical extinction and column density measurements 
for 7 supernova remnants (SNR), which yielded NH(cm~'^) 
= (2.22 ±0.14) X 10^^ Ay (mag). Later, using ROSAT obser- 
vations of 25 bright X-ray point sources as well as of 4 SNRs, 
iPredehl fc SchmittI (|l995l ) determined the shape and inten- 
sity of the soft X-ray halos around these sources, which they 
used in conjuction with dust halo models to measure the 
dust column density. In this study, they also found a relation 
between the X-ray derived hydrogen column density and op- 
tical extinction given by NH(cm~^) = (1.79 ± 0.03) x 10^^ 

Ay. 

The three analyses quoted above resulted in relations 
that disagree within their statistical uncertainties. This can 
be a result of (i) the low quality of early X-ray spectro- 
scopic data used; (ii) the systematic uncertainties between 
the methods employed, (Hi) formal uncertainties because of 
the small number of sources used, or (iv) strong variations 
in the ISM along the various lines of sight. In this paper, we 
take advantage of the most recent X-ray and optical obser- 
vations of SNRs with high quality instruments and compile 
the most complete set of extinction measurements towards 
these sources to redetermine the relation between the opti- 
cal extinction and the hydrogen column density. We choose 
SNRs because they do not sufi'er from significant intrinsic 
absorption that may contribute to the measured extinction. 
In addition, owing to the large number of observed sources 
and the high quality data from the latest generation of X- 
ray detectors, we will be able to address which of the above 
factors contribute to the variations in the Ay — Nh relation 
found in the previous studies. 

In Section 2 we summarize the methods used for the 
measurement of the hydrogen column density and optical 
extinction and present results towards 21 supernova rem- 
nants. In Section 3, we report the resulting Ay — Nnrelation 
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Table 1. Optical extinction and X-ray measured hydrogen column densities as collected from literature. 



SNR Name 


Nh 


Modeli 


Ref. 


E(B-V) 


Av 


Method 


Ref. 




(10^2 cm-2) 






mag. 


mag. 






SNR GO.0+0.0 


14.7+1.2 


TP 


(1) 


9.355+0.645 


29+2 


NS 


(2) 


SNR G004.5+06.8 


0.52 


SE 


(3) 


0.806+0.29 


2.5+0.9 


Fell Ratio 


(4) 


SNR G6.5-0.1 


0.47+0.1 


TP 


(5) 


1.15+0.15 


3.57+0.47 


Hc/H^ 


(6) 


SNR G13.3-1.3 


0.55+0.45 


TP 


(7) 


0.15 


0.47 


H^/H^ 


(7) 


SNR G53.6-2.2 


0.78+0.4 


TP 


(8) 


1.15+0.15 


3.57+0.47 


H<,/H^ 


(6) 


SNR G54.1+0.3 


1.6+0.1 


PL 


(9) 


2.581+0.226 


8.0+0.70 


NS 


(10) 


SNR G67.7+1.8 


0.515+0.195 


TP 


(11) 


1.7+0.3 


5.27+0.93 


Hc/H^ 


(12) 


SNR G69.0+2.7 


0.3+0.1 


PL 


(13) 


0.8 


2.48 


Hc/H^ 


(14) 


SNR G74.0-8.5 


0.023+0.01 


TP 


(15) 


0.08 


0.25 


NS 


(16) 


SNR G109. 1-1.0 


1.12+0.3 


AE 


(17) 


0.995+0.205 


3.15+0.65 


Hc/H^ 


(18) 


SNR Gill. 7-02.1 


1.25+0.03 


BB 


(19) 


1.613+0.129 


5.0+0.40 


SII ratio 


(20) 


SNR G116. 9+0.2 


0.79+0.12 


TP 


(21) 


0.871+0.161 


2.70+0.5 




(22) 


SNR G119. 5+10.2 


0.28+0.05 


TP+PL 


(23) 


0.410+0.132 


1.27+0.41 


Extinction map 


(24) 


SNR G120.1+1.4 


0.435+0.045 


TP 


(25) 


0.6+0.039 


1.86+0.12 


NS 


(26) 


SNR G130. 7+3.1 


0.416+0.08 


TP+PL 


(27) 


0.68 


2.11 


Hc/H^ 


(28) 


SNR G132.7+1.3 


0.43+0.25 


TP 


(5) 


0.710+0.040 


2.201+0.124 


Hc/H^ 


(29) 


SNR G184.6-5.8 


0.36+0.004 


PL 


(30) 


0.50+0.060 


1.55+0.186 


Ly« Absorption 


(31) 


SNR G260. 4-3.4 


0.454+0.049 


2BB 


(32) 


0.839 


2.60 


NS 


(33) 


SNR G263. 9-3.3 


0.0259+0.0001 


2BB 


(34) 


0.074 


0.38 


Hc/H^ 


(35) 


SNR G296. 5+10.0 


0.1+0.01 


2BB 


(36) 


0.161 


0.50 


Hc/H^ 


(37) 


SNR G327.6+14.6 


0.07+0.01 


SI 


(38) 


0.11 


0.34 


HI/GC 


(39) 


SNR G332. 4-00.4 


0.7+0.2 


2BB 


(40) 


1.516+0.29 


4.70+0.90 


Fell Ratio 


(4) 



^ Model abbreviations: TP: Thermal Plasma, PL: Power-Law, 
Synchrotron Emission, SI: Shocked ISM 



BB: Blackbody, AE: Absorption Edge modeling, NS: Nearby Stars, SE: 



Re ferences 
|Rho 



1 



)ISakano et al.l ||2004|): (2 ) |Pred"ehr 

Pctrj lll99?''- "'''^ 
(2008); (l l)lHui 

iKats uda ct al. (2008); (16) Parker (196 

J 20) Hurford & F cscn (19 96); (21) C raie et al.l (Il997l) 
I 



); (12 



Truernperl ll 19941 ); (3) [Reynolds et al.l (|2007|); f4)l01iva et al.l (Il989l); (5) 

(6) 'Long et al.' (Il99ll); (7)IS eward et al.' ("1995'): (8)'Sakcn ct al.' ('l995ll; f9 )lLu et al.l (1200 2^; (10) ' Koo et al] 
Becker (2008); (12) iMavrom atakis ct al. (2001); (13) Li ct al. (2005) ; (14)lHester fc Kul karni (1989|); (15) 

(17) iDurant fc van Kerkwiik (2006); (18) Fcscn & Hurford (1995"); (19) Hwang ct al. foW- 
(22) Fcscn ct al. (1997); (23) SJanc ct al. (1997); (24) Mavromatakis ct al 



2000); (25) Wa rren et a l .l l|2005): (26) 'Ruiz-Lapuentel_(l2004l) ; (27) iGotthel f ct al. (20o|); (28) lFesen eral] (Il988i) ; (29) Fcscn et al 

003) ; 



lUggS); (30) .Massaro et al.l mO^\ (31) ISoUerman et al.l ll200d); (3 2) iHui 



(120071 ); (35) [Manchester et al.l (Il978[) ; (36) lDe Luca et al.l (|2004|); (37) [Ruig (119831 ) 

( Il995l) ; (40) lDe Luca et al.l (l2006h . 



Becker (20061) ; (33) [Gorensteinl (fl975l) ; (34) .Manzali et al.l 
(38) lAcero et aO (|2007| ); (39) iRavmond et al.l 



and discuss any systematic uncertainties that may affect this 
result. 



2 EXTINCTION AND X-RAY ABSORPTION 
MEASUREMENTS 

Of the 24 3 supernova remnants i nclude d in the extensive 
catalog of iGuseinov et al.l (|2003l . l2004al lbh. approximately 
143 sources have been observed with the Chandra, XMM- 
Newton, or Suzaku satellites. The high quality spectral data 
obtained with these satellites allowed for a more precise mea- 
surement of the hydrogen column density by modeling both 
the X-ray continuum and the line features in the source 
spectra. We compiled the hydrogen column density mea- 
surements for these 143 sources (see Table 1 for the relevant 
references) . 

We then searched for independent measurements of 
the optical ext i nction towar ds all the SNRs listed in the 
iGuseinov et al. I (|2003l . l2004al lbh catalog. Unlike the case of 
X-ray extinction, our search revealed that there is a lack of 
reported measurements of optical extinction or reddening: 
Out of the 243 SNRs, we were able to obtain only 22 inde- 
pendent Ay measurements. Luckily, this set of 22 sources is 



a subset of the 143 sources, for which we have high-resolution 
measurements of the X-ray extinction with the new X-ray 
satellites. 

We discuss here the data as well as the methods that are 
used to determine the hydrogen column density and optical 
extinction. This will be important in quantifying systematic 
uncertainties in the correlation between optical and X-ray 
extinction. 



2.1 X-ray Hydrogen Column Density 
Measurements 

Hydrogen column density measurements of supernova rem- 
nants in the X-rays are performed by modeling their spec- 
tra, usually in the 0.2 — 8.0 keV range. Extinction is de- 
termined by comparing an intrinsic source spectrum, such 
as a blackbody, a power-law, or a thermal plasma emit- 
ting bremsstrahlung radiation, with the observed spectrum. 
The model parameters inferred from spectral fits are of- 
ten correlated with the hydrogen column density, especially 
when the analyzed data cover a small range in the X-ray 
b and and the energy resol ution is low (see, e.g.. Figure 2 
of lPredehl fc SchmittllTgQSI ). Moreover, when the spectra of 
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extended sources are observed with detectors of low angular 
resolution, the intrinsic variations of the source spectra that 
cannot be resolved bias the results. 

Observations with the Ch andra, XMM-N e wton , 
and Suzaku sat e llites (see, e.g.. iPaerels fc KahnI l2003l : 
lJuett et all 12004 lJuett et all 120061 ) vield significant im- 
provements in both of these areas. The increased angular 
resolution of the telescopes make it possible to model 
small regions within the SNR and detect local changes in 
the intrinsic spectral properties. At the same time, the 
increased energy resolution and signal to noise ratios lead to 
better modeling of both the continuum spectra and of line 
features. In the most favorable cases, e.g., SNR G109. 1 — 1.0, 
X-ray absorption edges of heavier elements, such as Mg and 
Ne, can be directly detected, completely eliminating the 
dependence of the column density on the assumed intrinsic 
X-ray spectrum of the source (iDurant fc van KerkwiikI 
I2OO6I ). All of these lead to better constrained hydrogen 
column densities compared to earlier studies. 

In the majority of the SNRs in our study, the X-ray 
spectra were fit with a thermal plasma model, where emis- 
sion lines were also taken into account whenever possible 
(see Table 1). The remaining spectra were modeled either 
by a power-law or a combination of two blackbodies. In or- 
der to make the analysis as free from the assumed continuum 
model as possible, we did not exclude from the study the re- 
ported hydrogen column density measurements where other 
continuum models were used. 

For the spatially resolved SNRs, hydrogen column den- 
sities measured towards different regions of the remnant can 
show variations that are larger than the statistical errors of 
each measurement. In these cases, we adopted the average 
of these values with an error that accounts for the observed 
scatter. Finally, in the cases where the central point sources 
can be resolved, we also used the Nnmeasurement obtained 
by modeling the emission from the central neutron star sur- 
face either by a blackbody, a power-law, or a hydrogen at- 
mosphere model. 



2.2 Optical Extinction Measurements 

Determining the optical extinction towards SNRs is more 
challenging than the hydrogen column density measure- 
ments in the X-rays and demands high signal to noise ratio 
spectra. Nearly all of the methods involve measuring the 
reddening using the intensity ratio of two emission lines and 
converting the reddening into an optical extinction. 

One of the well known and reliable methods of measur- 
ing the extinction is based on the Balmer decrement, which 
involves the intensity ratio of the observed Hq (6563 A) 
and (4861 A) emission lines. The observed relative in- 
tensity is compared to that expected for a gaseous nebula 
with typical temperature and electron density. Since this 
ratio depends very weakly on the physical conditions of the 
nebula, the theoretical ratio can be calculated with minimal 
uncertainty. Thus, the comparison allows a measurement of 
the r eddening, and hence, of the optical extinction (see, e.g., 
lOster brock 1989,,; Lequeux. ,2005 ). The majority of the opti- 
cal extinction measurements presented here are found using 
this method since both of the lines lie in the optical range 
and, in most cases, are strong enough to be resolved. 

One other frequently used method involves measur- 



ing the SII multiplet ratios (|Milleij [ 19681 ) in the infrared 
(~10320 A) and blue (~4068A). As with the Balmer decre- 
ment method, this ratio also depends only very weakly on 
the temperature and density in the remnant. However, the 
disadvantage of this method is the fact that it is not always 
practical to perform spectral observations in both of these 
wavelength regions with high enough spectral resolution. 

In a third method, the two most prominent IR transi- 
tions of Fe[II] (1.6435 and 1.2567/im) that arise from the 
same upper level can be used. Their intensity ratio does not 
depend strongly o n the temperature and density of the gas 
(|01iva et al.lll989l ). allowing for a determination of the ex- 
tinction. In our sample, the optical extinction towards SNR 
G332. 4-00.4 and SNR G004.5-I-06. 8 has been measured with 
this method by l|Qliva et al.l 19891). using a theoret ical ratio 
of 1.36 for these lines (Nussbaumer fc Storevl[l98i) . 

A final method involves using nearby stars with known 
distances to estimate the extinction towards the SNR. In Ta- 
ble 1, we denoted the four optical extinction measurements 
obtained with this method as "nearby stars" . 

In Table [T] we present the compilation of the optical 
extinction data. For the cases where only the reddening 
E(B— V) measurement was reported, we converted these to 
an extinc tion (Ay) assuming a Ay = 3.1E{B — V ) relation 
(see, e.g.. lFitzpatric^l2004l : [Savage fc Mathidl 19791 ). Finally, 
in the absence of any reported systematic or statistical er- 
rors, we assigned a 15% uncertainty on the measurements 
for the purposes of the fit only, which is similar to the av- 
erage errors given for the other measurements. These cases 
are left without errors in Table 1. 



3 RESULTS AND DISCUSSION 

We plot in Figure 1 the optical extinction, Av, for 22 su- 
pernova remnants against their hydrogen column density, 
measured from their high energy resolution X-ray spectra. 
Figure 1 includes a typical 15% error on the optical extinc- 
tion for those remnants that did not have reported errors on 
this measurement (see Table 1). We performed a linear fit 
between these quantities and obtained the best-fit relation 

NH(cm~^) = (2.21 ± 0.09) x 10^^ Ay, (1) 

where the errors correspond to 1—a statistical uncer- 
tainty. We present the best-fit line in Figure 1. This rela- 
tion can also be expressed as NH(cm~^) — (6.86 ± 0.27) x 
10^^ E{B — V){mag) between the reddening and the hydro- 
gen column density. 

We then investigated whether the continuum models 
that are used to fit the X-ray observations to measure the hy- 
drogen column density affect the results found for the Ay — 
Nnrelation. Specifically, we looked into the power-law con- 
tinuum models, which are known to give higher values for 
the hydrogen column density than others. To this end, we 
compared the Ay— Nh relation that we found from fitting 
the entire sample to that derived when the three power- 
law continuum fits in the sample were excluded. We found 
that the resulting relations do not differ significantly from 
each other either in the slope or in the errors. Specifically, 
excluding the subset of observations where the X-ray contin- 
uum was modeled with a power-law results in the relation 
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Nh= (2.29 ± 0.11) X 10^^ X Av, while all the data yield 
Nh= (2.21 ± 0.09) X 10^^ X Av. 

Another possible systematic uncertainty is introduced 
by the use of solar abundances when measuring the hy- 
drogen column density from X-ray spectra. It has been 
shown that different assumptions on the elemental abun- 
dances (e.g., solar vs. ISM) in the interstellar medium can 
lead to ~ 5% var iation in the hydrogen column density 
I Wilms et al.ll2000l ). The Nnvalues used in our study could 
be subject to similar systematic uncertainties. 

Comparing our result with previous studies, we find 
th at it is in very go od agreement with the relation given 
bv lGorensteinI l| 19751 ) ^n{cm-'^) = (2.22 ±0.14) x 10^^ Av, 
while it differs at the 2>a level from the result obtained by 
iPredehl fc SchmittI (|l995l ) l<ln{cm-'^) = (1.79 ±0.03) x 10^^ 
Av. This difference can be attributed to a number of causes: 
ROSAT's narrow bandpass, which does not allow tight con- 
straints on the intrinsic spectral shape; the authors' assump- 
tion that the intrinsic spectra of all their sources are power 
laws; possible intrinsic absorption; and the uncertain prop- 
erties of the optical counterparts, which affect the assumed 
extinction values. 

Despite the improvements on the number and quality of 
the observations, there may still be a number of systematic 
errors and selection effects present in our study. A larger, 
homogeneous sample of supernova remnants, studied with 
both the last generation X-ray satellites as well as with sen- 
sitive optical spectrographs will help refine the studies of 
the interstellar medium and the relation between the optical 
extinction and the hydrogen column density in the Galaxy, 
hence the dust to gas ratio. Increasing the number of sources 
in the sample may also allow for a search for variations in 
the relation towards different lines of sight. 
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